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RAPID CONVERGENCE TO FREQUENCY 

FOR 
SUBSTITUTION TILINGS OF THE PLANE 

JOSE ALISTE-PRIETO, DANIEL CORONEL AND JEAN- MARC GAMBAUDO 



Abstract. This paper concerns self-similar tilings of the Euclidean plane. 
We consider the number of occurrences of a given tile in any domain bounded 
by a Jordan curve. For a large class of self-similar tilings, including many 
well-known examples, we give estimates of the oscillation of this number of 
occurrences around its average frequency times the total number of tiles in the 



domain, which depend only on the Jordan curve. 



1. Introduction 
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& Quasicrystals are alloys that show long range order but possess symmetries that 

prevent them from being crystals. From their spectacular experimental realization 
in the early 80's |SBGC84] to their very recent discovery as natural objects in the 
Kamtchatka mountains [BSYL09J, quasicrystals have been the subject of very active 
research, whose application and interest extend far beyond the scope of solid state 
physics. The first examples observed were rapidly quenched alloys of Aluminum 
t-H and Manganese exhibiting icosahedral symmetry. In these quasicrystals, atoms are 

known to appear with a given frequency. This means that for a large ball B of 
radius R, the ratio of the number N s of atoms inside B corresponding to a specific 
t ■ atomic element s (for example Aluminum in the first known cases) to the total 

C^ number N of atoms in the ball B tends to a limit v s as R tends to infinity: 

, ~~ l ,. N s 

I' hm — = v s . 

For such examples, the following questions arise naturally: 
!-h (1) to estimate the oscillation of the ratio N s /N around its limit v s , that is 

to say, to give upper bounds for the speed of convergence to of the ratio 
\N S /N — u s \ as R goes to infinity; 
(2) to extend these estimates to domains whose closure are not Euclidean balls. 
We now make these questions more concrete by using aperiodic tilings of the Eu- 
clidean space to model quasicrystals (see |Bcl03]), for details and precise statements 
of the results, see Section [2j Consider a tiling T of <i-dimensional Euclidean space 
R d made with isometric copies of a finite set of tiles {pi, . . . ,p n }, where the p^s 
are homeomorphic to a closed ball in M. d . Let 7bea (large) set in M. d also homc- 
omorphic to a closed ball and let dV be its boundary. For each i e {1, .. . ,ri}, 
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we let N{V,T,Pi) denote the number of isometric copies of pi in V, N{V,T) the 
number of tiles of T, and C(dV, T) the number of tiles of T that intersect dV. 
If V is a ball, then it is well-known that there are tilings such that for each i in 
{1, ...,n}, the ratio N(V,TiPi)/N(V,T) tends to a well-defined limit Vi as the 
radius of V tends to infinity (see for instance GS87, Sol97, LP03]). In fact, this is 
usually a consequence of the unique ergodicity of the associated dynamical system 
(see for instance [LMS02] ). However, in general, very little is known concerning an 
upper bound for the quantity \J\f(V,T,Pi) — v i .N{V 1 T)\ 1 specially when V is not 
a Euclidean ball (see |LP03[ lAPCj for the case when V equals a ball or box) . In 
particular, for an aperiodic tiling, can we hope to obtain an estimate as strong as 
the one we can formulate for periodic tilings which reads: 

(*) MV,T, Pi ) - *VV(V,T)| < KC(dV,T) 

for some K > 0? 

In this paper, we answer this question in the affirmative for a large class of self- 
similar tilings in two dimensions. Self-similar tilings, which are probably the most 
studied examples of aperiodic tilings (see for instance [Rad94, Sol97, AP98, Rob04| 
ICS06| and references therein, examples are given in Section [5]), are associated with 
a (substitution) primitive square matrix M having non-negative integer elements. 



By Perron- Frobenius theory (see for instance |Sen81[|HJ90] , see also Subsection 3.3 
for more details), we know that there exists a largest positive real eigenvalue /x, the 
remainder of the spectrum being in a ball centered at with radius smaller than 
ti. We denote by r(M) the modulus of the second largest eigenvalue of M, that is, 

r(M) = max{|ry| | r\ ^ fi is an eigenvalue of A/}. 

The simple condition r(M) < ^/Jl ensures that the above estimate (*) holds true for 
the associated self-similar tiling. If r{M) = y/JL and if all eigenvalues with modulus 
r(M) are semi-simple (their algebraic and geometric multiplicities coincide, see for 
instance [IJ90J for more details), then the estimate (*) has to be relaxed to: 

(**) |JV( V, T, Pi) - UiN{V, T) I < KC(dV, T) ]nC{dV, T) 

for some K > 0. Some explicit examples of self-similar 2-dimensional tilings satis- 
fying all these conditions will be given at the end of the paper. 

These results are more explicitly stated in the next section. Ideas for the proofs 
have been inspired by J. Peyriere's 1986 paper |Pey86| and can be related to B. 
Adamcezswki's work |Ada04j . 

2. Definitions and Main result 

2.1. Definitions. Let A be a closed subset of the Euclidean plane K 2 . A tiling of 
A is a countable collection T = (tj)j&j of closed subsets of A that cover A and have 
pairwise disjoint interiors. The sets tj are called tiles and, in this paper, all tiles 
are supposed to be homeomorphic to the closed unit ball in K 2 , see Section [5] for 
examples. 

Let £ be a group of isometries on the plane that contains the group of transla- 
tions. A tile q is £ -equivalent to a tile p (or is a £-copy of p) if q is the image of p 
by an isometry in £. Let V = {pi, . . . ,p n } be a finite collection of tiles. A tiling T 
of A C R 2 is £ -generated by V if every tile in T is ^-equivalent to some tile in V '. 
The set of all tilings of A that are £ -generated by V is denoted by ilg^-p(A). When 
A = R 2 , we write fls,v instead of Qg ^(IR 2 ). 



Speed of convergence for self-similar tilings 3 

Let A > 1. Given a subset U C R 2 , let XU := {Ax | x G U}. A substitution rule 
(with dilation factor A) is a collection S = (Si)" =1 , where 5^ belongs to ftevi^Pi) 
for all i S {1, . . . , n}. Let iS be a substitution rule. On one hand, the dilation by A 
induces a natural map from f2£ t -p(A) to ^^(AA), where XP = {\p \ p E V}. On 
the other hand, S induces a map from fls t \-p(AA) to r^.-p (AA), which is defined 
by subdividing the tiles of a tiling in fig A-p(AA) according to the substitution rule. 
For details see Subsection |3.2| The composition of these two maps yields a map 
1-s.A '■ ^f.-p(A) — > £lf -p(AA). If A = M 2 , then 1$ := 1-sm 2 is a self-map of Qe,v 
and is referred to as the substitution map. A tiling T is called self-similar if it is a 
periodic point for 1$. A tiling T in ^ls,v is admissible for S if it belongs to 

k>0 

The substitution matrix associated with S is the n-by-n matrix Ms — (mi ,-)* j, 
where for each !,je{l,...,n}, the coefficient m it j is the number of £ -copies of pi 
in Sj. Recall that a matrix M is primitive if there exists n > such that all the 
elements of M n are positive. The substitution rule S is primitive if its substitution 
matrix is primitive. In this paper, all the substitution rules considered will be 
primitive. 

2.2. Main result. Let V = {pi, . . . ,p n } be a finite set of tiles, S a substitution rule 
with dilation factor A > 1 and T an admissible tiling for S. Given a Jordan curve 
T in R 2 bounding a topological closed disk A, we let £(T, T) denote the number 
of tiles of T that intersect T, Af(A,T) be the number of tiles of T included in A, 
and for each ie {1, . . . , n}, we define Af(A, T,Pi) to be the number of £-copies of 
Pi contained in A. 

The following theorem, which constitutes the main result of this paper, provides 
estimates on the oscillation of the number of occurrences J\f(A 7 T,Pi) around an 
average frequency, which only depends on the Jordan curve T bounding A. 

Theorem 2.1. Let V = {p±, . . . ,p n } be a finite collection of tiles, S a substitution 
rule with dilation factor A > 1 and primitive substitution matrix Ms- There ex- 
ist positive numbers V\ , . . . , V n depending only on the substitution matrix Ms that 
satisfy: 

(i) If r(Ms) < A, then there exists K > such that, for every Jordan curve T 
bounding a topological closed disk A and every tiling T in Qs, *"e have: 

\Af(A,T, Pi ) - Uitf(A,T)\ < KC{T,T). 

(ii) Ifr(Ms) = A and all the eigenvalues with modulus A are semi-simple (that 
is, their algebraic and geometric multiplicities coincide), then there exists 
K > such that, for every Jordan curve T bounding a topological closed 
disk A and every tiling T in Qs, we have: 

\N{A,T,Pi) - Vl U{A,T)\ < K£(T,T)\n£(T,T)). 

Remark 2.2. Notice that Theorem \2.1\ is given for a quite general setting: 

(1) On one hand, we do not require the tiles in V to be polygons, they need not 
even be disks with piece-wise smooth boundaries. In fact, the boundary of 
the tiles could be Jordan curves with infinite length and positive Hausdorff 
dimension, as is the case for Rauzy tilings that we shall discuss later. 
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(2) On the other hand, we do not require the standard finite pattern condition 
used in tiling theory (see [KPOOjj. This allows us to deal with self-similar 
tilings which have "fault lines" along which tiles can slide past one another. 
This was the case for the self-affine tilings studied by N. Priebe Frank and 
L. Sadun [PFS09] . 



Remark 2.3. For a given substitution, the estimates given by Theorem \2.1\ may 
depend on the group £ being considered. This is the case of Penrose tilings, see 
Section^ for details. 



Remark 2.4. Theorem 2.1 can be easily extended to colored tilings (for an example, 
see the square tilings in Section [50. 

3. The tools 

3.1. Jordan curves and locally finite tilings. In this section, we provide some 
simple combinatorial estimates for the number of tiles of a tiling that intersect a 
Jordan curve. First we must define some notation. Denote by B x (r) the closed ball 
of radius r around x in M 2 . Given a subset U CM. 2 , define 

fmin(U) := sup{r > | there exists x £ U s.t. B x (r) C U}. 

To a tiling T of M 2 , we associate 

r T = inf{ r min (t) | t € 7} 

and 

R T = sup{diam(i) | t e T}/2. 

The following lemma shows that if rj- is positive and Rj- is finite, then there is 
a uniform bound on the number of tiles intersected by balls of a prescribed radius. 

Lemma 3.1. Suppose that T is a tiling with Rj- < +oo and rj- > 0. Then every 
ball of radius 2Rj- intersects at most Kj- := \16Rj-r^- J tiles of T, where [—\ 
stands for the integer part. 

Proof. Fix i£l 2 . It is clear that each tile of T intersecting B x (2R-j-) is included 
in B x (4R-j-). Since each one of these tiles contains a ball of radius r-j~, comparing 
area yields r^-nN < IQttR^-, where N is the number of tiles intersecting B x {2Rf) 
and the conclusion follows. □ 

By virtue of the previous lemma, we say that a tiling T is locally finite if Rj- < 
+oo and rj- > 0. The next result provides an estimate for the diameter of a 
simple curve in terms of the number of tiles of a locally finite tiling that the curve 
intersects. 

To simplify notation, it is convenient to identify curves with their images. Let 
7 : [0, 1] — > M 2 be a simple curve. The identification induces an order on the image: 
Let y,y' in 7 with y = j(s) and y' = j(t). If 7 is open, then y < y' if and only if 
0<s<t<l;if7is closed, then y < y' if and only if < $ < t < 1. 

Lemma 3.2. Let T be a locally finite tiling and 7 be a simple curve. Then, 

diam( 7 ) <2i? r £( 7 ,T), 

where £(7, 7") denotes the number of tiles that 7 intersects. 
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Figure 1. Proof of Lemma 3.2 Construction of the sequence 
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Proof. By compactness, there exist x, x' G 7 with x < x' and d(x, x') = diam(7). 
We construct a sequence of points (xj)f_ as follows. Fix e > and set xo = x. To 
construct X\, if a:' belongs to B Xo (2Rj- + e), then set £ = 1 and 27 = x'. If x' does 
not belong to B Xg (2Rj- + e), then define #1 by 

Xi := max{z e 7 | z e B Xo (2Rj- + e)}. 

It is clear that d(x\,Xo) = 2Rq--\-s. This construction can be extended inductively 
to obtain a sequence x = xq < x\ < . . . < xg-\ < xi — x' in 7 such that for all 

ks{i,...,i-iy. 

• d(x k -i,x k ) = 2R T + e, 

• {z£j\x k <z<x'}n U^} B Xm {2R r + e) = {x k }, 

and such that d{xi-i,xt) < 2Rj- + e (see Figure [TJ. 

For each k G {0, ... ,1 — 1}, choose a tile tj-(fc) in T that contains x^. Since 
d(xk,Xk>) > 2Rj- + e for all fc, fc' in {0, ...,£ — 1}, it follows that £,•(&) is not 
the same tile as tj(k') unless k = fc'. Hence, the number C("f,T) of tiles in T 
that intersect 7 is at least £ (see Figure [I]). Moreover, it is easy to check that 
£(2Rj- + e) > d(x, x') = diam(7). Thus, using the estimate for £(7, T) we get 

diam( 7 ) < (2i? r + e)£(7,r), 

and the conclusion can be obtained by letting e go to 0. □ 

Given a Jordan curve Y and two tilings T and T', the next lemma compares 
the number of tiles of T that Y intersects with the number of tiles of T' that Y 
intersects. 

Lemma 3.3. Let Y be a Jordan curve and T be a locally finite tiling. Then, for 
every locally finite tiling T' such that Rj-i > Rj- and Kj-i < C(Y,T') < C(Y,T), 
we have: 

C(Y,T') < (2K T > + 1)^C(T,T), 



where Kf is the constant defined in Lemma\3.1 
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Proof. For every point y G M 2 , we define 

C y :={t£T'\tnB 2R , T (y)^<b} and C y :=U teCy t. 

We construct a collection y = {yi}? =1 (with p to be determined) of points in T such 
that the collection {C yi }? =1 covers T. Fix any point of T as y\. Now suppose that 
2/1, . . . , yj have been chosen such that, for alH £ {1, . . . , j} and k £ {1, . . . ,i — 1}, 
the point y t does not belong to C Vk . There are two cases to consider. Either the sets 
{C yi } J i=0 cover T, in which case we set p = j and the construction is completed; or 
they do not cover T, in which case we choose any point of r\U^ =0 Ut e c . £ to be j/j+i 
and continue iterating the construction. Observe that £(r,T') is finite because T' 
is locally finite and T is compact, and in each iteration, we add at least one tile of 
T' that intersects T to the area covered by {C Vi }l_ . Hence, the construction stops 
in finitely many steps. 



On one hand, by Lemma 3.1 each C Vi contains at most Kq-i tiles of T'. It 
follows that 

(1) C(T,r)<pK T >, 

and since K-f < C(T,T'), we deduce that p > 2. 

On the other hand, fix i £ {1, . . . ,p}. Hence, there exists j £ {1, . . . ,p} such 
that yj, which belongs to T, does not belong to B Vi (R-j-' — Rt)- Suppose that 
Hi < Uj (the other case is analogous), then there is a point y[ > j/i in T with 
d{yi,y[) = Rt> - Rt sucn that the arc [y t ,y'i\r = {z £ V \ y t < z < y-}, which 
joins yi and y[, is contained by the ball B y .(R-f> — Rt)- Since i was arbitrary, using 



Lemma 3.2 yields 

Rr> -Rt< diam([^,^] r ) < 2R T C([y i ,y' i \ r ,T) 

for all i £ {1, . . . ,p}. Combining all these inequalities, we obtain 

v 

(2) p(r t > - Rt) < 2R T J2^([yuy'i\r,r)- 

From the construction of y, it is clear that d(yi,yj) > 2Rf for all i < j, and, in 
particular, the collection B = {B yi {R-f — Rt)W=\ is pairwise disjoint. Moreover, 
the distance between two different balls in B is greater than 2Ry- It follows that no 
tile of T may intersect more than one ball in B. Hence, no tile of T may intersect 
more than one arc [j/i, 2/i]r- Thus, combining Q and (IT]) we get 

(3) c(r,T) < 2K T , Rr £(r,r). 

txrfi — ttf 

To finish the proof, fix c > arbitrarily and consider the following two cases. First, 
suppose that R T ,/R T < 1 + c. Since £(I\T) < C(T,T), it follows that 

(4) c(T,T')<(l + c)^C(T,T). 

Rt> 

Now suppose that Rf /Rt > 1 + c. It is not difficult to check that 

{\ + c){R T , -R T ) > cR T >- 
Replacing this inequality in p| we get 

(5) C(T,T')<2K T ^^l + ^jC(T,T). 
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Combining ([5| and Q yields 

An easy computation shows that the last bound is optimal when c = 2Kq-i and the 
conclusion follows. □ 

3.2. Hierarchical sequences. In this subsection, we recall the concept of "hirar- 
chical sequence" of a tiling. This notion have been used in the context of self-similar 
tilings satifying the standard finite pattern condition by many authors, see for in- 
stance KPOCk Here, we extend its use to admissible tilings. 

Let S be a substitution rule with dilation factor A. For each I € N, let \ l V = 
{X l p | p G V} and consider the map Di : &£,\i-i-p — > ^s,\'v defined by Di{T) — 
XT '. Also consider the decomposition maps Si : Q £ \ l v ~ * ^s.x'- 1 v defined by 

Si(T) = {giX^Si) \g€£,ie{l,..., n},g(X l Pl ) e T} 

for each I e N* . By definition, I s — S\o D x . 

It is not difficult to check that Is is onto when restricted to Sis- This implies 
that for each admissible tiling T, there is a sequence (T l )i £ n of tilings, called a 
hierarchical sequence of T, such that T° — T and T /_1 = Si(T l ) for all I <G N, 
that is, each tile of T l can be decomposed into tiles of T /_1 according to the 
substitution rule \ l ~ 1 S. This sequence is constructed as follows: Set T° = T and, 
for each I > 0, inductively choose T l € 1$ 1 (T /_1 ) and then set T l — \ l T l . It is not 
difficult to check that X l - l S 1 (X i - l T') = Si(V) for every V in tt £ ^ v . It follows 
that T'" 1 = Si{T l ) for all I e N. 

Remark 3.4. Notice that all tilings in Qs are locally finite and that for each tiling 
T and each I > 0, T l is also locally finite and satisfies: 

r-j-i = A rf, R-fi = A Rf, and thus K-p = Kj-. 

Proposition 3.5. LetT be an admissible tiling for S. For every topological closed 
disk A there is a finite collection Ao, . . . , A TO _i of closed subsets of A such that: 

(*) u|z^ _1 Aj = Ujgj A *j where J A = {j \ t 3 C A}. 
(ii) All A/ 's have pairwise disjoint interiors. 
(Hi) For each I € {0, . . . , m — 1}, A; is a union of tiles in T l , which does not 

contain a tile in T l+1 ■ 
(iv) A does not contain a tile of T m ■ 

The collection A , . . . , A m _! is called a hierarchical decomposition of the closed disk 
A. Moreover, ifT is the Jordan curve bounding A, then 

(6) Y n -'- 1 <^ J C(r,7*) and N{A h T l )< ||Af|| 1 £(r,T i+1 ) 

rr 

for all I € {0, . . . , m — 1}, where \\M \\\ is the maximum absolute column sum of the 
substitution rule M . 

Proof. Choose m to be the smallest integer such that no tile in T m is included in 
A and set 

Pm-i = {* € T™- 1 | t C A} and A m ^ 1 =U tePm _ 1 t. 
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Applying the recursion 

Pi-i={teT 1 - 1 \tC A\A,} and Aj_i = A; U (U te pk-i*). 

for I G {1, . . . ,m— 1}, we obtain a sequence of sets Pq, . . . , P m _i and A , . . . , A m _ 1; 
and it is straight forward to check that the sequence A; = U te p,i for all I £ 
{0, . . . ,m — 1} satisfies properties (i) to (iv). 

To check the first inequality in mty, observe that, on one hand, Lemma 3.2 implies 



diam(A) < 2R-fiC(T,T l ). On the other hand, since A contains a tile of T" 
it follows that 1rx m < diam(A). The inequality is now obtained by replacing 
R-j-i = \ l Rj- and r-j-m-i = X m ~ l rf in the previous inequalities. 

It remains to check the second inequality in (pi). From the construction of the 
Pi's, it is clear that each tile in Pi is contained in a tile of T l+1 that meets V (which 
is the boundary of A). Since each tile of J~ l+1 is subdivided in tiles of T l according 
to the substitution rule, it follows that 

Af(A u T l )< ||M|| 1 £(r,T m ). 

□ 

3.3. Perron-Frobenius Theory. We now recall the basic Perron-Frobenius the- 
ory that we will need in the sequel. For proofs see JHJ901 Chapter 8] . Let M = Ms 
be the substitution matrix associated with S. Since S is primitive, the matrix M 
is a primitive matrix with non-negative integer coefficients. The Perron-Frobenius 
Theorem states that the largest real eigenvalue p > of M , which is called the 
Perron eigenvalue, is simple and greater than one. Moreover, there exists a right 
eigenvector v T = (vi, . . . , v n ) and a left eigenvector w = (wi, . . . , w n ) such that v 
and w have positive coefficients and (v,w) — 1. We denote by r{M) the modulus 
of the second largest eigenvalue of M, that is, 

r(M) = max{|?7| | r\ ^ /i is an eigenvalue of M}. 

Recall that an eigenvalue is called semi-simple if its algebraic multiplicity is 
equal to its geometric multiplicity. The following proposition, or more precisely the 
corollary below, will be crucial for the proof of our main result. The first part is a 
well-known consequence of the Perron-Frobenius Theorem, see for instance [HJ9Q1 
Theorem 8.5.1]. The second part is less well-known, but can be easily deduced 
from the first part by using the Jordan Canonical Form of M after recalling that 
the Jordan blocks associated with semi-simple eigenvalues are diagonal. For / > 
and i,j £ {1, . . . ,n}, m\ denotes the (i, j)-element of the matrix M l . 

Proposition 3.6. For every p > r(M) there exists K — K(p) > such that 

(7) Imlj-ViijtwjlKKft, 

for all I > 0. Moreover, if the eigenvalues of modulus r(M) are semi-simple, then 
there exists K > such that ufh holds with p — r(M). 



Corollary 

there exists 


3.7. Let v t = 
a constant K 


V i/J2j V j f° r al1 ' 

> such that for 


£ 

all 


i € 


{!,' 


n}. 


Then, 
n}, 


for 


all p 


> 


r(M), 


(8) 




n 

\rn\j — Vi 2_^ m 
fc=i 


k,j 


< 


Kp l 
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c) I = m 



Figure 2. Construction of the hierarchical sequence for a given 
A. The tiling in gray T° is a Penrose tiling (see Section pi). The 
tiling in black is T l , the dark-gray region is A; and the light-gray 
one is A/ + i. In this example, m = 5. 



for all I > 0. Moreover, if all the eigenvalues with modulus r(M) are semi-simple, 
then there exists a constant K > such that ([8]) holds with p = r{M). 

Proof. An easy computation shows 



J2 m kj ^ \ m h -Vil^Wj 



k=l 



E, 



2J \v k H l Wj 
fe=l 



%j\ 



The conclusion now follows from applying Proposition |3.6| twice. 
4. Proof of the main result 



□ 



Let V = {pi, . . . ,p n } be a finite set of tiles and S a primitive substitution rule 
with dilation factor A > 1. Suppose that T is an admissible tiling for S and T is a 
Jordan curve in K 2 bounding a closed disk A. 



10 



J. Alistc-Pricto, D. Coroncl, J.-M. Gambaudo 



The idea of the proof is as follows. Fix i <G {1, . . . , n} and let Vi be defined as in 
Corollary 3.7 First, it is not difficult to check that the number off-copies of tj after 
applying the substitution I times to ti is exactly m\ ■ . It follows that the density of 
the tile ti is Vi. Second, we consider a hierarchical sequence (T l = {t l j)j>o)i>o of T 
and the hierarchical decomposition {Aq, . . . , A m _i} of A constructed in Proposition 



3.5 We use the decomposition of A to estimate A/"(A,T,Pi) and jV(A,T)- Since 



the sets A; have disjoint interiors, we have 



(9) 

and 

(10) 



Af(A,T, P i) = Y,tf(Ai,T,Pi) = J2 E -^(4,r, Pl ) 



1=0 



'=0 tJCAi 



771—1 

1=0 tic A, 



Next, multiplying (10 1 by Vi and then subtracting (pi), we get 

m— 1 

(11) M(A,T, Pl ) - ^(M(A,T)) = Y, E M(ti>T,Pi) - v^^T). 

1=0 t'CAi 



Suppose that ft is a £-copy of A'p,t for some k € {1, 
mkandJV(4T)=E < 



, n 



}. Then7V(4,T,ft) = 



fc . Therefore, applying Corollary 3.7 to (11), we get 
m — 1 ?7i— 1 



|AA(A,T,^)-^AA(A,T)| < E E K P l -- E^^' 7 "')^' 

i=0 ('CA, 



/=() 



where if > 0, and either A > p > r(M) if r(M) < A or p = A if r(Af ) = A and the 
eigenvalues of modulus r(M) are semi-simple. Applying Proposition 3.5 to the last 
inequality, we obtain: 



(12) 



MA,T, Pi ) - vMA,T)\ < K\\M\U J2 C(T,T l+1 )p l 



1=0 



We want to apply Lemma 3.3 to give an upper bound of C(T,T l ) in terms of 
£(r,T) for all I € {0, . . . ,m}. Since for each I > 0, the tiles of T l+1 are tiled by 
tiles of T l , it follows that £(r,T 7 ™- 1 ) < . .. < ^(^T 1 ) < £(I\T). It is easy to 
check that K -p — Kj- for all I > 0. 

Next, we define l € {0, . . . ,m} as follows. If £(T, T" 1 " 1 ) > K T , then ^ := m. 
Else, we define ?o to be the minimal Z € {0, . . . , m — 1} such that £(r, T'°) < -K"r- 
Observe that in both cases C(T, T la ~ x ) > K T for all I £ {0 . . . , l - 1}. 



We split the sum in ( 12 ) into two parts 

Jn-l 



m— 1 



J 



(is) ^£(r,rV + E £ ( r > r V 

1=0 l=l 

and deal with each part separately. For the first sum, applying Lemma |3.3[ we get 



R T 



C{T,T)<{2K rl +l)-^C{T,T) 
tt-j-i 
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for all I € {0, . . . , Iq — 1}. Since Rf/Rr 1 — A l f° r every I > 0, it follows that 



Jo-l 



Jo-l 



(14) 



^£(r,rV < (2ifr + iK(r,r)X; (f) 



;=o 



z=o 



To estimate the second sum in (13), we suppose that Iq < m (otherwise, the sum 
is zero). Since p < X and £(T,T l ) is decreasing, we have 



m— 1 



m — 1 



£ £(r,rV < £(r,r<°) £ a* = c(t,t Lo ) 



i=i 



i=h 



A-l 



From Lemma ^ and the fact that C(T,Ti ) < K T , we get A"'" 1 "' < K T R T /r T . 
Denote N T := (XK r R r /r r - 1)/(A - 1). It follows that 



m — 1 



(15) 



^ £(r,rV < n t x 1o c{t,t 1 °) 



l=la 



If /o > 0, then from Lemma 3.3 and (15), we get 

m— 1 



(16) 



J2 C{T,T l )p l < (2K T + l)N T C(T,T). 



i=i 



and (16), we get 



It is clear that (16) also holds when ^o = 0, since (2Kj- + 1) > 1. Combining (14) 



n -i 



(17) M*,T, Pi ) - «/*AT(A,T)| < (2K T + 1)£(I\T) ( £ (j ) + A> | . 



,i=0 



Finally, we deal with the different cases for p. In the first case, we have r(Ms) < 
p < X and it follows from (17) that 



\Af(A,T, Pi ) - PiAf(A,T)\ < (2K T + 1) f-^ + Nr) £(T,T), 
which finishes part (i) of the main result, since Rf, Tj- and Kj- are constant in fig. 



In the second case, we have p = X and then from ( 17 1 and the estimation of m > £q 



given in Proposition 3.5 we obtain 



\Af(A,T, Pi ) - ViN{A,T)\ < (2K T + l)C(T,T) (log A C(T, T) + N T ) , 
where Nj- = Nj- + 1 + log x (R-r/r-r)- Thus, there exists K > such that 



\M{A,T,Pi) -ViN{A,T)\< KC{T,T)hxC{T,T), 

as soon as T meets at least 2 tiles. This completes the proof of part (ii) of the main 
result. 
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5. Examples 

5.1. Penrose tilings. Penrose tilings are among the most known examples of self- 
similar tilings, see for instance [Sen95, GS87, Pcn80l lA"P98] . Here, we follow the con- 
struction of Anderson and Putnam [AP98] . Consider two isosceles triangles t\ and 
£2 in R 2 , where the vertices of t\ have coordinates (sin(7r/10), 0), (— sin(7r/10), 0) 
and (0, cos(7r/10)) and those of £2 have coordinates (sin(37r/10), 0), (— sin(37r/10), 0) 
and (0, cos(37r/10)). Both triangles are equipped with decorations (arrows) on their 
edges as shown in Figure [31 

Let £ denote the group of translations in M 2 and Pp nroso be the set of 40 
triangles obtained by rotation of one of the two triangles t\ or £2 with its prescribed 
decoration by an angle kir/10, where k is in {0, 1, . . . , n — 1}, and let rip cnrosG be the 
set of tilings of M 2 made with translated copies of triangles in "Ppcmosc such that their 
tiles meet full-edge to full-edge and the decorations on overlapping edges coincide. 
Elements of f2p cnroso are called Penrose tilings. They can also be constructed by 
using the substitution <Sp en rose described in Figure [3] The dilation factor of <Sp cnrosc 
is the golden mean A = (1 + \/5)/2. 

The substitution matrix Mp onrosc is a 40 x 40 non negative primitive matrix 
(see |AP98j for details) . It is well-known that the Perron eigenvalue of Afp cnroso is 
equal to A 2 = (3 + a/5)/2, and that A and —A are the only eigenvalues of modulus 
r(Mp onrose ). Straight-forward computations show that both of them are semi-simple 



with multiplicity 3. Thus, we can apply Theorem 2.1 (part (ii)), and obtain: 

Each triangle pi in Ppenroso has a well-defined frequency i>\ > 0, and there exists 
K > such that, for every Jordan curve T and every tiling T in r2p cnrose , we have: 

\tf(A,T,Pi) - ^AA(A,T)| < KC(T,T)HjC(T,T), 

where A is the closed disk bounded by T. 

Now, Let £ be the group of all isometries (direct and indirect) of M 2 . Up to these 
isometries, there are now only 2 types of Penrose tiles (the fat and thin triangles). 
The associated substitution matrix reads 

2 1 
1 1 

and (3 + V5)/2 and (3 — v5)/2 are its eigenvalues. Thus, we can apply Theorem 



2.1 (part (i)) to obtain: 

\M{A,T, Pi ) - vMA,T)\ < KC(T,T). 



5.2. The square and table tilings. Squares and table tilings were studied by 
Robinson [Rob04, Rob99J. First, we consider the square tilings. Let 'P S q Uare be 
the set of four unit squares respectively decorated with symbols p,q,r and s. The 
square substitution rule 5 square is described in Figure 4] The dilation factor is 
A = 2. In this case, the substitution matrix Af squarc reads 

(2 1 l\ 

12 10 

12 1 

\1 1 2) 
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Figure 3. Local and substitution rules to construct Penrose tilings 



s_P 

q p 







r s 

»- 

r q 



q -*- 


q q 

P r 


H- 


P r 
s s 



Figure 4. Square table substitution 



Figure 5. The table substitution 



The eigenvalues of M square are A 2 , A and 0. It is easy to check that A is semi-simple 



with multiplicity 2. Applying Theorem 2.1 (part (ii)) (extended to deal with the 
decorated case), it follows that each square t £ Square has a well-defined frequency 
v t > and there exists K > such that for any Jordan curve T bounding a close 
disk A and any tiling T in fls , we have: 

<J o '-'square > 

W(A,T,t) - v t Af(A,T)\ < KC(T,T)\n(C(T,T)). 

Now we consider the table tilings. Let "Ptabie = {Pi,P2} where pi is a vertical 
domino and p2 is a horizontal domino. The substitution rule <St a t>ie ^ s described in 
Figure [5j In this case, the dilation factor A is 2 and the associated substitution 
matrix reads 

(2 2^ 

Its eigenvalues are A 2 and 0. Applying Theorem 



2.1 



(part (i)), it follows that each 
domino pt has a well-defined frequency Vi > 0, and there exists K > such that, 
for any Jordan curve Y bounding a close disk A and any tiling T in ^s tablc , we 
have: 

\Af(A,T, Pl ) - ViN{A,T)\ < K C(r,T) . 
The square and table tilings can be identified (see |Rob04] for details) by iden- 
tifying tiles as described in Figure [6j 
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HE 



Figure 6. Identification of table and square tilings 





Figure 7. The pinwheel substitution 



This identification can be used to count tiles in the different tilings. Consider for 
instance the unit square marked with a p. In every tiling in Qg s uarc , the squares 
marked with p and r always appear as in Figure [6j Thus, for any closed Jordan 
curve r bounding a close disk A we have: jV( A, T, p) — Af( A, T, p\ ) is smaller than 
£(r,T), where p\ is the horizontal domino in 'Ptabic- Hence, using the estimates 
for the table substitution gives better estimates than the estimates for the square 
substitution. It is natural to wonder why? The answer simply relies on the fact 
that the square substitution system is less efficient than the table one in the sense 
that the estimates of the error for the nth substitution of the square p, which is a 
square of size 2™ , yields an error term of size A™ = 2™ . This error term is actually 
going to be canceled by terms corresponding to the other squares of size 2™. In our 
computation, the error term for a collection of squares of size 2 n is bounded from 
above by the sum of the error terms for each square of size 2" and thus we cannot 
see these cancellations which actually occur. 



5.3. The pinwheel tilings. Pinwheel tilings were introduced by Conway and 
Radin Rad94 . Let £ be the group I + of direct isometries of ]R 2 , and Ppinwheei 
be the set of two right triangles which are isometric (but with inverse orientation). 
Each triangle has a hypotenuse oflength 1 and sharp angle a satisfying tana = 1/2. 
The substitution rule Spinwheel is described in Figure [7] and has a dilation factor of 

'2 3\ 
„ „ ) It has two eigenvalues A 2 = 5 



The substitution matrix is the 2x2 matrix 



and — 1. Applying Theorem 2.1 it follows that each triangle ti in "Ppmwhoci has a 
well-defined frequency Vi > 0, and there exists K > such that, for any Jordan 
curve r bounding a close disk A and any tiling T in Q$ lnw h..i : 

MA,T, Pi ) - v t M{\T)\ < KjC(T,T). 



5.4. Rauzy tilings. Rauzy tilings (see |Rau82j for details) are tilings made from 
a set 'Prauzy of three (topological) disks ri,V2 and r%, where r*i is a disk whose 
boundary is a Jordan curve (with positive Hausdorff dimension), r 2 = Ari and 
r 2 = A 2 r 3 where A ~ 1.84 is the unique real root of the polynomial X 3 — X 2 — X — 1 
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(see |Fog02| ). The substitution matrix is the 3x3 matrix 




The largest eigenvalue of the substitution matrix is A 2 . All other eigenvalues have 
modulus smaller than 1. Applying Theorem |2.1| it follows that each Rauzy tile pi 
has a well-defined frequency i/* and there exists K > such that for every Jordan 
curve r bounding a close disk A and every tiling T in £1$, we have: 

\N(A,T, Pl ) - Vitf(A,T)\ < KL{T,T). 
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